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Рассматривается задача о критических нагрузках сжатой пря- 
моугольной пластины, содержащей непрерывно распределен- 
ные источники собственных напряжений. Анализ задачи про- 
водится на основе модификации системы нелинейных уравне- 
ний Кармана для больших прогибов упругих пластин с дисло- 
кациями и дисклинациями с различными вариантами краевых 
условий. Введением замены для функции напряжений задача 
сводится к исследованию двух задач: линейной краевой зада- 
чи относительно функции напряжений, вызванных внутрен- 
ними источниками и системы нелинейных уравнений относи- 
тельно прогиба и функции напряжений, вызванных внешними 
сжимающими нагрузками, которая имеет тривиальное реше- 
ние. Классическая критическая нагрузка определяется как 
наименьшее собственное значение линейной краевой задачи, 
полученной линеаризацией системы нелинейных уравнений 
относительно тривиального решения. Рассматриваются четы- 
ре типа краевых условий: все края подвижно защемлены; все 
края шарнирно оперты; два противоположных края свободны 
от напряжений, а два других подвижно защемлены или шар- 
нирно оперты. Равномерно распределенные сжимающие 
нагрузки одинаковы на противоположных краях. Установле- 
но, что если мера несовместности является нечетной по одной 
переменной и четной или нечетной по другой переменной, то 
напряжения, вызванные только внутренними источниками, не 
приводят к потере устойчивости плоского равновесного со- 
стояния и не влияют на критические значения сжимающих 
нагрузок. 


Ключевые слова: упругая пластина, дислокации и дисклина- 
ции, критическая нагрузка. 
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Введение. В работе Л. М. Зубова [1] на основе уравнений Кармана построена система нелинейных уравнений 
равновесия упругих пластин, содержащих внутренние источники напряжений, вызванных дислокациями и дисклина- 
циями. Рассмотрена задача об изгибе мембраны вследствие релаксации вызванных дефектами внутренних напряже- 
ний, которая сведена к уравнению Монжа-Ампера. Найдено несколько точных решений о форме поверхности мем- 
браны, содержащей распределенные дисклинации. В [2] рассмотрена задача о влиянии внутренних напряжений на 
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прогибы и напряженное состояние круглой упругой пластины, испытывающей поперечное давление. Установлено, 
что наличие дисклинаций приводит к нелинейному увеличению прогиба. Исследованы задачи об устойчивости 
и послекритическом поведении нагруженной контурным давлением круглой пластины с непрерывно распределенны- 
ми дисклинациями. Найдены вызванные дисклинациями осесимметричные закритические формы изгиба пластинки, 
которые существуют при отсутствии внешней нагрузки. Установлено, что при переходе пластинки из плоского состо- 
яния в изогнутую форму уменьшается величина внутренних напряжений. 

В работе В. А. Треногина [3] был развит метод Ляпунова—Шмидта в операторной форме для нелинейных уравне- 
ний в банаховых пространствах. В работе [4] Л. С. Срубщик и В. А. Треногин исследовали задачу о влиянии малой 
поперечной нагрузки на выпучивание и послекритическое поведение пластины произвольной формы под действием 
параллельных осям координат сжимающих краевых усилий. 

В работе Рейсснера [5] рассмотрена задача о влиянии нелинейно-упругого основания на выпучивание и началь- 
ное послекритическое поведение безмоментного плоско-напряженного состояния в случае бесконечной пластины с 
малыми геометрическими несовершенствами, а в случае тонкой сжатой пластины строго выпуклой формы со свобод- 
ным краем при дополнительном действии малой поперечной нагрузки эта задача решена Л. С. Срубщиком [6]. В рабо- 
тах [7, 8] исследована задача о влиянии малой поперечной нагрузки на послекритическое поведение прямоугольной 
гибкой пластины, лежащей на нелинейно-упругом основании и равномерно сжатой в продольном направлении, а в [9] 
эта же задача рассмотрена для пластины с дислокациями и дисклинациями. С помощью операторного метода Ляпуно- 
ва-Шмидта определено количество решений, соответствующих новым формам равновесий, и для каждого из них 
строятся асимптотические представления. 

1. Постановка задач 

1. 1. Система уравнений равновесия и краевые условия. Пусть тонкая упругая прямоугольная пластина, со- 
держащая в плоском состоянии поля непрерывно распределенных краевых дислокаций и клиновых дисклинаций, 
находится под действием малой нормальной нагрузки интенсивности ЕС(Х ‚У ) и внешнего краевого усилия, состав- 


ляющие которого вдоль осей Х и У равны соответственно Р и О. Тогда система уравнений равновесия записывается в 
виде [1, 10] 


РА? = [+,"|+ЕС 
А + ЕЙ, /2 = Ей (1.1) 
р = ЕВ? /12( - у?), [Фи] = Ух Фу — Му ху + МУХ 
Здесь 
ш(Х,У) — скалярная мера несовместности, которая выражается через плотности дислокаций И дисклинаций, 


У’(хХ и ) — прогиб пластины, $(х ‚У ) — функция напряжений Эри, й — толщина пластины, у — коэффициент Пуас- 


сона, Е — модуль Юнга, Х, У — прямоугольные координаты с началом в центре пластины и осями, параллельными ее 
краям, & — малый числовой параметр. 


Систему (1.1) будем исследовать с краевыми условиями одного из следующих типов [9, 10] 

























































































У =Их =Ч ху =0, 4, =-Р/|Х|=а/2,|У| < 5/2 #2 
=, = Фу =0,4х =-0,У|=Ь/2,|Х| < а/2 
М =Ух =Ч ху = 0,4, =-Р/Х|=а/2,У| <Ь/2 аа 
М =Уу = Фху = 0,4 хх =-О/У| = Б/2,|Х| < а/2 
У =Ух = Фу = 0,4, =-Р/Х|=а/2|У| < 5/2 
учу х = Иру + (2 УМху = Чу = хх =0У| =Ь/2,|Х| <а/2 —. 
У =Мх = Ч ху =0, 4, =-Р|Х| =а/2,У| < 5/2 а 
Ух = Му + (2) хху = Чу = Фх =0У| =Ь/2,|Х| < а/2 
Краевые условия (1.2) отвечают свободному защемлению, а (1.3) — подвижной шарнирной опоре краев. 


Сжимающие усилия равномерно распределены по краям Х =+а/2 и У=+/2 соответственно. Условия для и =р/2 


в (1.4) и (1.5) соответствуют свободному краю. В этих случаях сжимающие усилия приложены только к краям 
Х ==а/2. 
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Пустъ Ч=Ф-Р у" й 2-0хХ ь / 2. Тогда функция Ф удовлетворяет однородным краевым условиям 
Фуу =Фу, =0 при || =а/2, Фуу =Фух =0 при У =Ь/2. Заменяя эти условия на эквивалентные Ф =Фу =0 при 


|| =а/2, Ф=Ф, =0 при |" =Ь/2 , перепишем систему уравнений (1.1) и краевые условия (1.2) — (1.5) в виде 


Е: +0; = Ф.М +66 а 


А?Ф + ЕНИ/, У /2 = Ей 








(1.7) 





М =И, =Ф=Ф, =0,Х|= а/2,|У| < 5/2 
\ =И =Ф=Ф, =0,|У|=5/2.|Х| <а/2 











(1.8) 





М =Ихх =Ф=Фх =0,Х|=а/2,У|<Ь/2 
\ =, =Ф=Ф, =0,У|=Ь/2,|Х|<а/2 











(1.9) 








М =Их =Ф=Фх =0,|Х|=а/2,У|<5/2 
ИУ хх = Му + (2-х =Ф=Ф, =0|У| =Ь/2,|Х| <а/2 








(1.10) 








И =И/х =Ф=Фу =0,Х|=а/2,У|<Ь/2 
|. ЧУ х = Му + (2-у) ду =Ф=Ф, =0У|=Ь/2|Х|<а/2 
В случае краевых условий (1.9) или (1.10) считаем, что О=0 в (1.6). 
1.2. Переход к безразмерным переменным. В (1.6) — (1.10) перейдем к безразмерным переменным по формулам 
Х =ах,У = Бу М(Х,У)= их, у) /у,Ф(Х ,У) = РЕ(х, у) 
Р=Бр/Ь?,О = Ра (ва? 6 (х,У)= паз (х,у)/(ь*] (1.10) 
а=Ь? /а? и = ру? /п? 
С учетом (1.11) система (1.6) преобразуется к виду 
Аи + рам +а\, = о, Е [9 
о $ <12, |У<1/2, (1.12) 
А. ( )=()» +9), Аа =А.Аа 
Введем числовой параметр г такой, что (ох, У) = Не У), где в (х, у) фиксированная функция. Обозначим 
через Е решение краевой задачи 
№ Е=(х,У), (1.13) 
Е=Е, =0,|4 =2, |< 2 а 
Е=Ё, =0,=12х <У2 
Тогда решением задачи АЕ = Вх, У) с краевыми условиями (1.14) будет функция Ё, =7Ё\. Назовем пере- 
менную Г параметром интенсивности напряжений, вызванных внутренними источниками пластины. 


Положим Е = 7Н + Е. Учитывая (1.11) — (1.14), выводим краевые задачи 


А+ рам» +4, = ам, Е]+ паи, Е +88 15) 
АЕ, /2=0, |4 <1/2, |< 12, | 
=, =0, = 1/2, || </2 
и, 02 не 
ш=и, =0, || = 2, <2 


= =0,|х| =12, 1/2 
| вх 0.4 в 


и=и =0,| =1/2.х <1/2 


ый ба 2 
| и, =0,=/2,]5|<У (1.18) 


м, мм = ии, + (2-ммм, =0у|=И2х < 1/2 
У У» хх 
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&=1/2,| | <2 


=. =0, 
"У + Уби’,„ = Му + (2 = Ум у = 0, | = у 2. 





(1.19) 





«<1/2 

Функция Е в (1.15) должна удовлетворять краевым условиям (1.14). 

При & =0 краевая задача (1.15) с каждым из краевых условий (1.16) — (1.19) при всех значениях параметров 
р, 4, гимеет тривиальное решение (и», Е,)= (0,0). 

2. Критические нагрузки выпучивания основного равновесия 

2.1. Исследование линеаризованных задач. В работе [9] с использованием результатов [11-13] показано, что 


критические нагрузки выпучивания основного равновесия определяются из линейной краевой задачи на собственные 
значения 






































Ао-+а®,, —т[,Е]= — Ра (2.1) 
с параметрами 4, г и одним из следующих краевых условий 
=®,=0,х=12, 1/2 
Ф=О, х=12, |< В 
Ф=®, =0, У =1/2, || <1/2 
= =0,|х| =1/2, 1/2 
Ф=® =0,|х|=1/2,|у| < В 
®=ю =0,|у|=1/2| <1/2 
=®,=0,1х| = 12, 1/2 
Ф=®, х=12, |<! 2 
Уи =, +@-удю и, =0,|]=И2,М <1/2 
=®„ =0,х| = 2, 1/2 
@=ю =0, |4 =1/2,|у| < (2.5) 
0) + У; =® у + 2- У) =0, У = 2, х 12 
Здесь в случае краевых условий (2.4) или (2.5) в уравнении (2.1) полагается а=0. 
Введем линейные операторы А и В 
Ао = А ®+4®,, _ го, Е |, Во = —0© (2.6) 


Задачи на собственные значения (2.1) с краевыми условиями (2.2—2.5) в случае г=0 исследованы, где пока- 
зано, что оператор А формально самосопряжен и ограничен снизу [14]. Формальная самосопряженность оператора А 
в случае г=*0 следует из равенства 


Го, Рау = По орам 


где функция & удовлетворяют одному из краевых условий (2.2)-(2.5), а функция Е— условиям (1.14). Ограничен- 
ность снизу оператора А в случае г 0 следует из неравенства 


1 о, Е рахау Поло 


которое следует из неравенств Фридрихса [14]. 
Так как в задаче на собственные значения (2.1) оператор А формально самосопряжен и ограничен снизу, а 


<с 











[оо 4 
о 








оператор В положительно определен [13], то можно применять результаты $839, 40, 44 из [15] и свести задачу к ми- 
нимизации функционала 





Ко)= т ы ‚ (м,ъ)= аа р (2.7) 


где функция © удовлетворяет одному из краевых условий (2.2)-(2.5). Приходим к следующему алгоритму [14, 16]. 
Пусть известны первые собственные значения р <р.<..р,| и соответствующие собственные функции 
61, 62,....0„_1 задачи (2.4) при заданных значениях параметров 4 и г. Задаем начальное приближение 0) к соб- 
ственной функции @®,. Для Т= 0,1... выполняем следующие действия 
1) вычисляем (. + 1) — е приближение к собственному значению р», 
рб") = 50) (2.8) 


2) находим 5 из уравнения 
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Аб = р во (2.9) 


п 


при соответствующих краевых условиях из (2.5)—(2.8); 


3) ортогонализируем решение 5) к первым п-1 собственным функциям по энергии оператора В 
б : п-1 й 
(9) = 56) У (85, (Во о) (2.10) 
р 


4) вводим функцию о) 


о) в 5/5 





п ) 








|6] = (ах [вх У) (2.11) 


<=, где & — за- 





в 
Этот процесс повторяется до тех пор, пока не будет выполнено условие сходимости | ре и р® 


данная погрешность. 
Критическое значение параметра х› определяется как наименьшее собственное значение задачи 
2 — 
Доо- рас „а, = по, Е | (2.12) 
при заданных значениях параметров р и 4 с каждым из краевых условий (2.22.5) (в случаях (2.4), (2.5) полагается 
4=0). Если в (2.12) р=а=0, то при г=х пластина теряет устойчивость плоского равновесного состояния под 


действием только внутренних источников напряжений. 
Для решения задачи (2.12) применяется минимизация функционала (2.7), где операторы А и В определяются по 
формулам 


Ав = А 0+ раю +аюу, Во =а®,Е | (2.13) 


у? 
Здесь оператор А ограничен снизу [15], а знакоопределенность оператора В зависит от свойств функции ЕЁ, т.е. от 
свойств меры несовместности 0% У). 

В силу симметрии краевых условий (2.2)—(2.5) относительно начала координат при условии четности или не- 
четности функции Ё по каждой переменной собственные функции @ задач (2.1) и (2.12) являются соответственно 
четными или нечетными по каждой переменной. Отсюда следует, что если функция Е является нечетной по одной 


переменной и четной или нечетной по другой переменной, то для оператора В в (2.13) справедливо 
(Ве, о) = а] ©, Е рахау = а в, Е ахау = 0 (2.14) 
о О 


и для оператора А в (2.6) имеем 
(Ас, ©) = (о +а®„- по, Е 1 ®)= (о И о] (2.15) 
Из (2.14), (2.15) следует, что задача (2.12) имеет только тривиальное решение © =0 ‚ав задаче (2.1) собственные 
значения не зависят от параметра г. Это означает, что если мера несовместности и (х У) является нечетной по одной 


переменной и четной или нечетной по другой переменной, то напряжения, вызванные внутренними источниками, не 
приводят к потере устойчивости плоского равновесного состояния и не влияют на критические значения сжимающих 


нагрузок р и 4. 
Если функция Е такова, что оператор В в (2.13) положительно определен, то задача (2.12) разрешима методом 
минимизации функционала (2.7), т.е. критическое значение г = к) существует. 

В данной работе для каждой конкретной функции Ё|, являющейся решением задачи (1.13), (1.14) при задан- 
ной мере несовместности и (х У) существование решения задачи (2.12) определялось численными экспериментами, 
т.е. проверялась сходимость алгоритма (2.8)—(2.11), где вместо параметра р, используется г,, а операторы А и В 
определены формулами (2.13). Полученные при заданных значениях параметров р = ру, а = 4% собственное значение 
то и собственная функция © задачи (2.12) сверялись с собственным значением р» и собственной функцией ©» зада- 
чи (2.1) при заданных значениях г=х, 4=40. Во всех случаях сходимости алгоритма (2.8)—(2.11) наблюдалось 
р» = Ро и ©, = ®‹ с погрешностью менее 0,001. 


2.2. Результаты численных расчетов. Для краевых условий (2.3) решения уравнения (2.1) при г=0 выписы- 
ваются явно [10] 


Рт.п = Гдббий +и?} ваз? | Двдзи?) (2.16) 


Вестник Донского государственного технического университета 2016, №1(84), 43-51 





„и = зи (лт(х + 0,5) (пи(у +0,5)) (2.17) 


где т, п — натуральные числа. Для остальных краевых условий при г=0 и для всех краевых условий при г 0 зада- 
ча (2.1) решалась по алгоритму (2.8)—(2.11) с применением разностного метода и квадратурной формулы Симпсона. 
Формулы (2.16), (2.17) использовались в тестовых расчетах при отладке программ. Из (2.16) следует, что при 


г=0 иа=л” ( — 40? у а,0<а<!/2 минимальному собственному значению риши = Ри! = Р›1 = п? (ба+2) 0 отвечают 
две собственные функции &, | — четная по обеим переменным, и ®›,— нечетная по переменной х, четная по пере- 


менной у. 


Рассмотрим случай постоянной меры несовместности, т.е. в задаче (1.13) и (х, у)= 1. Тогда решение задачи 
(1.13), (1.14) — функция Е является четной по обеим переменным, что подтверждается численным решением этой 
задачи разностным методом. На рис.1 для случая квадратной пластины приведены графики функции Ё| и её вторых 
производных Ни, Ну, Ву › Которые позволяют судить о распределении напряжений, обусловленных внутренними 


источниками. 





Е Е хх 


Рис.1. Графики функции Ё! и её вторых производных Ёу» Еуу» Еьху 


С помощью численных экспериментов найдены критические значения параметров нагрузок р, 4, г которым 
отвечают одна или две собственные функции при различных значениях геометрического параметра а=Ь? |. а? и крае- 
вых условиях (2.2)(2.5). Отметим основные результаты для прямоугольных пластин при а =1/4, 1, 4. Для случая 
г=к, р=О, 4=0 (пластина теряет устойчивость плоского равновесного состояния под действием только внутрен- 
них источников напряжений) установлено, что для закрепленной по краям платины (условия (2.5), (2.6)) к, >0,а для 
пластины с двумя свободными краями (условия (2.4), (2.5)) к) < 0. На рис.2 приведены графики первой и второй соб- 


ственных функций задачи (2.12) с условиями (2.4) для квадратной пластины при р=0, 4=0. 





Рис.2. Графики первой и второй собственных функций задачи (2.12) с условиями (2.4) 
для квадратной пластины 


Пешхоев И.М. О критических нагрузках сжатой упругой прямоугольной пластины 





В табл.1 приведены значения первого и второго собственных значений и, г› задачи (2.12) с условиями (2.4) 
для квадратной пластины при различных значениях параметра р (а = 0). 


Таблица 1 
Собственные значения п, г› задачи (2.12) с условиями (2.4) 





Р "1 72 

0 -3794 | 3719 
7,73 |-3282 | 3334 
15,45 | -2741 | 2910 
23,17 | -2123 | 2424 
30,89 | -1371 | 1812 
38,62 | 0,11 | 590 



































В табл. 2 приведены критические значения и для р=4=0 при различных значения параметра @ =Ь?/ а, 


а =1 и краевых условиях (2.2) - (2.5). 


Таблица 2 
Критические значения % для р=а=0 




















Краевые условия | а=1/4 | а=1 | а=4 
(2.2) 40807 | 4354 | 2550 
(2.3) 27856 | 2595 | 1740 
(2.4) -16997 | -3794 | -3088 
(2.5) -10166 | -2120 | -1708 

















Заключение. 

1) Для рассматриваемого в работе случая симметрии краевых условий относительно начала координат уста- 
новлено, что если мера несовместности является нечетной по одной переменной и четной или нечетной по другой 
переменной, то напряжения, вызванные внутренними источниками, не приводят к потере устойчивости плоского рав- 
новесного состояния и не влияют на критические значения сжимающих нагрузок. 

2) Для постоянной меры несовместности на примере квадратной пластины при различных значениях сжима- 
ющей нагрузки численными экспериментами показано существование положительных и отрицательных критических 
значений параметра интенсивности напряжений, вызванных внутренними источниками (табл.1). 


Автор благодарит профессоров Л. М. Зубова и Л. С. Срубщика за помощь в работе. 
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